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We investigate the quasiparticle interference in the heavy Fermion superconductor CeCoIns as 
direct method to confirm the d-wave gap symmetry. The ambiguity between d^y and dj,2_y2 sym- 
metry remaining from earlier specific heat and thermal transport investigations has been resolved in 
favor of the latter by the observation of a spin resonance that can occur only in d^2_y2 symmetry. 
However these methods are all indirect and depend considerably on theoretical interpretation. Here 
we propose that quasiparticle interference (QPl) spectroscopy by scanning tunneling microscopy 
(STM) can give a direct fingerprint of the superconducting gap in real space which may lead to a 
definite conclusion on its symmetry for CeCoIns and related 115 compounds. The QPI pattern for 
both magnetic and nonmagnetic impurities is calculated for the possible d-wave symmetries and 
characteristic differences are found that may be identified by STM method. 

PACS numbers: 74.25.Jb, 71.27.-|-a,72.15.Qm 



I. INTRODUCTION 

In recent years the spectroscopic imaging scanning tun- 
neling microscopy (SI-STM) has become a powerful ex- 
perimental tool for studying the local electronic proper- 
ties of various superconductors^ It is well known that 
the Fourier transform of STM (FT-STM) data or quasi- 
particle interference (QPI), can be used to elucidate the 
nature of the many-body states in novel superconduc- 
tors, in particular those having quasi-two-dimensional 
electronic structure. In the presence of impurities, elastic 
scattering mixes two quasiparticle eigenstates with mo- 
menta ki and k2 on a contour of constant energy. The 
resulting interference at wave vector q = k2 — ki reveals 
a modulation of the local density of states (LDOS). The 
interference pattern in momentum space can be visual- 
ized by means of the SI-STM^ 5. 

In layered cuprates the analysis of the QPI has pro- 
vided details of the band structure, the nature of the su- 
perconducting gap, or other competing orders- ^. It was 
shown that a magnetic-field dependence in quasiparticle 
scattering interference patterns is sensitive to the sign of 
the anisotropic gap^iJill-, Recently, such QPI effects have 
also been studied in the iron based superconductorJi^, to 
find the order-parameter symmetrjli^HlSI ^ case where 
we have multi-band superconductivity with possible sign 
reversal of the order parameter between the electron and 
hole pockets. 

Applying this technique in the heavy Fermion sys- 
tems and considering QPI in these materials is one of 
great interest, in particular because there are numer- 
ous unconventional superconductors with competition of 
magnetism and superconductivity as well as non-Fermi 
liquid behaviour. One of the most difficult issues in 
the heavy Fermion superconductors like the 115 com- 
pounds CeMIns (M = Co, Ir, Rh) is the identification of 
the symmetry of superconducting order parameter. In 
the present context of QPI theory the microscopic origin 
of the unconventional superconductivity is not an issue. 



For CeCoIns various proposals based on spin fluctuation 
theorjff^I^ and composite pairin^l^l have been advanced. 

Commonly several candidates for the gap function are 
proposed and their compatibility with temperature and 
field dependence of thermodynamic and tr anspo rt quan- 
tities is used to discriminate between thenpQISU jj^ 
respect it is instructive to recall the previous discussions 
on the gap symmetry in this compound. Firstly the Pauli 
limiting behaviour of the upper critical field-^ and ob- 
served KnighlP^ shift proves the spin singlet nature of 
the gap suggesting d-wave pairing. Originally field-angle 
resolved thermal conductivity^'^ and specific healP^ ex- 
periments which probe the node structure gave conflict- 
ing results of dx2-y2 and dxy gap symmetries respectively. 
Then the observation of a pronounced spin resonance^ 
at a;r/2Ao = 0.65 in the superconducting state with in- 
elastic neutron scattering (INS) gave strong evidence for 
the former. Namely detailed calculations of the spin re- 
sponse with realistic Fermi surfacc^^ show that the res- 
onance can appear only for the d^^-y^ symmetry but 
not in the dxy case. Further field-angle resolved specific 
heat measurements at even lower temperatur^^ finally 
also concluded on dx2_y2 gap symmetry. These inter- 
pretations however all depend considerably on theoret- 
ical model features and approximations. For example, 
the quasiparticle relaxation rate in the vortex state for 
the transport properties influences the results of field- 
angle resolved specific heat measurements. In addition, 
the results of the INS experiments were also interpreted 
in terms of spin wave excitations which sharpen in the 
superconducting state due to effect of the gap on the nor- 
mal state Landau damping.l^Sl It would be preferable to 
have a more direct method that can provide a fingerprint 
of superconducting gap symmetry. In this work we show 
that QPI can indeed serve this purpose in heavy Fermion 
superconductor CeCoIns and other related 115 systems 
where the degree of the three-dimensionality is more sub- 
stantial than in cuprates or iron-based superconductors. 
In Sec. [IT] we present the theoretical framework of this 
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FIG. 1. (Color online) a) Calculated Fermi surface for 
CeCoIns using the band structure parameters defined in 
Ref. [Slj, b) first panel indicates a cut through the Fermi 
surface for normal, and second and third panel indicate a cut 
through the Fermi surface for d^y and d^2_y2 gap symmetry 
respectively (solid red lines) for bias voltage lj — O.lAo. (at 
— first row, at fcz = O.Stt second row, and at fc^ = tt last 
row). Dashed (green) lines indicate the node- lines, ± denotes 
the sign of the superconducting gap, and are the typical 
scattering vectors defining the QPI pattern. 



Finally Sec. |IV|gives a brief sum- ^^^^^ hybridized quasiparticle {a±Ma) bands 



model. Then Sec. |III| presents the numerical results for 
the QPI and how to interpret them in terms of Fermi 
surface and nodal gap properties of the proposed order 
parameter candidates 
mary and outlook. 



II. THEORETICAL MODEL 

As a starting point we need a model for the elec- 
tronic structure which captures the essence of the heavy 



quasiparticle bands and their Fermi surface but is tech- 
nically still manageable for a T-matrix calculation of 
QPI in the superconducting state. The electronic struc- 
ture of CeMIns has been investigated by using tight 
binding models with hybrid isation for f-electrons and p- 
(conduction) electrond^nM]^ -j-j^g Ce-4f electron states are 
split by a large spin orbit-coupling (Ago ^ 0.4 meV^Hl) 
into upper i—7/2 and lower j=5/2 multiplets. One may 
therefore restrict to the lower one as done in Ref. [30] 
which is further split into three crystalline electric field 
(CEF) Kramers doublet states. Because the CEF split- 
ting energy is about three times the heavy quasiparticle 
band width {W ~ 4 meV) one may further restrict to the 
lowest CEF doublelP^ which has an effective pseudo-spin 
1/2. Then the Anderson lattice model Hamiltonian for 
the two hybridized (c,f ) conduction and localized orbitals 
which are doubly spin degenerate is given by 



^ = E ^kcLck. + + (cL/ka + h.C. 

kcr 

+T.Ufffltf^tfl'ih'i- (1) 



kk' 



where cj^.^ creates an electron with spin a in the conduc- 
tion orbital with wave vector k = {kxjkyjkz). Further- 
more, and ej|^ are effective tight binding dispersions 
of the conduction band and the renormalized dispersion 
for the / band respectively. /^^ creates the f electron 
with momentum k and pseudo spin cr, and Uff is its 
on-site Coulomb repulsion. Finally Vk is the hybridiza- 
tion energy between the lowest 4f doublet and conduc- 
tion bands which contains implicitly the effect of spin 
orbit and CEF term. A use the Anderson lattice model 
of Ref. [3T] is preferable here because it provides a conve- 
nient way to realistically model the heavy quasiparticle 
bands. For complementarity we mention that there is 
also an alternative method to investigate SI-STM, pro- 
posed in Refs. and |33J which starts from the Kondo 
lattice model where charge fluctuations are eliminated 
and 4f electrons are considered as fully localized. 

It is known thal!^ in the limit J7// — >■ oo where double 
occupation of the f-states are excluded, by defining the 
auxiliary boson one can find the mean field the Hamilto- 
nian as 



IMF — E 
±,k(T 



T~l-MF — 



± t 



1 



4± 



(2) 



with = V-^{l—nf) denoting the effective hybridisation 
obtained by projecting out double occupancies. Due to 
1 — n/^lVkis strongly reduced with respect to the sin- 
gle particle 14- Using the parameters defined in Ref. [31] 
for the above quasiparticle band structure we plot the 
corresponding Fermi surface (FS) in Fig. ([T a). The FS 
crossing originates only from one of the bands, E^, while 
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remains well above the Fermi energy. Therefore, in 
the following analysis of the QPI it does not play a sig- 
nificant role. 



To supplement the normal state Hamiltonian we add 
the superconducting pairing term given by 

-Hsc - E (4,kT4.-H + ^■^) (3) 

k 

where is the superconducting gap function. 

The scattering of the quasiparticle by perturbations in 
the sample such as non-magnetic or magnetic impurities 
is responsible for the QPI which is believed to be mea- 
sured in SI-STM . We perform the analysis of such pro- 
cesses based on a T-matrix descriptiorP. In particular, 
we introduce an impurity term in the Hamiltonian 

'Himp — ^ ^ ("^kk''^o'(T' + •^ctit'S • CTaa') .kcr'^i .kV i 

kk'CTCT' 

(4) 

where J^^., and Jo-o-' represent the non-magnetic and the 
magnetic point-like scattering between the electrons re- 
spectively. Note that in this picture Ce does no longer 
carry a local magnetic moment because below the co- 
herence temperature its f-electron is incorporated in the 
itinerant quasiparticle states. The magnetic impurities 
may be other 4f ions with stable moment. Depending 
on the crystalline electric field splitting their coupling to 
quasiparticle states can become strongly anisotropic, e.g. 
of Ising type in the simplest case which will be consid- 
ered in the following. This means we set the quantization 
axis of the magnetic impurity spin along the z-direction 
and consider only the Sz spin component. Here, S refers 
again to the pseudo spin of the lowest Kramers doublet of 
the 4f impurity moment. At this stage we do not include 
scattering of quasiparticles by nearly critical collective 
spin fluctuations which may be important in CeCoIus 
because of its closeness to an antiferromagnetic quantum 
critical pointP^. 

Deflning the new Nambu spinor as = 

(a^ j^^, a+__k^, ^..j., a-^-ki), the Hamiltonian can 
be written as 

^ = II v^I/SkV-k + ^k^kk' v-k' (5) 

k kk' 

By introducing J^^., = 7; and JzzSz ~ 7', the above 
matrices /3k and U^k' are defined as 

/3k- {^^^)^{E+<Jz + Ai,<j,) 

+ (^'^y^)<»(^k^^+^kff-); (6) 

C^kk' = To (8) (7'ao + 70-2). (7) 

Here tXi are the Pauli matrices acting in spin space, Ti 
are the Pauli matrices in the orbital space, and Ti (Si (Ji 
denotes a direct product of the matrices operating on the 
4-dimensional Nambu space. 



FIG. 2. (Color online) The Quasiparticle interference due 
to non-magnetic impurity impurity scattering, for individual 
layers at k^c — 0.0, 0.27r, 0.47r, O.Gtt, O.Stt, tt, for normal state 
(first panel), superconducting state with dxy gap symmetry 
(second panel) and d^2_y2 gap symmetry (third panel). (For 
small bias voltage lu — O.OlAo and using Aq/W = 0.086 
where W is the total quasiparticle band width.). 
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Therefore in terms of the Nambu spinor the Green's 
function (GF) matrix in Matsubara representation, is ob- 
tained via Gkk'(''') — ~(^V'k(''')V'k'(0))' whence 

Guk'K) = G^K)[<5kk' + Tkk'K)G^,(c^„)], (8) 

where G^(lli„) — ^«w„ — /3k^ is the unperturbed prop- 
agator (bare GF) of the conduction electrons. Its poles 
are given by the cf-hybridized (due to Vk) quasiparticle 
excitations gapped by Ak which are contained in the 
/3k matrix. As a result of the hybridisation gap of order 
Vk the QPI interference pattern will depend strongly on 
the energy or bias voltage within the gap range. 
Solving the Dyson equation for the T-matrix 

Tkk' {y>n) = C^kk' + ^ t^kk"Gk„ (w„)Tk"k' (t^ri)^ (9) 
k" 

the LDOS can be obtained from imaginary part of the full 
GF, via analytic continuation iujn — > a; + iO^ according 
to 



N^{uj) = — Im Tr 

TT 



( 2 ^ 



^N^{lo) + 6N^{u). 



(10) 



The Fourier transform of the fluctuations, 6N^{ijj), can 
be obtain using Eq. dSl) in following form. 



1_ 



-Im Tr 



(11) 



E 

k 



To 



- ) Gk (W„ )Tkk+q (t^n ) Gk+q ( W„ ) 



This quantity is called QPI and qualitatively is propor- 
tional to the convolution of density of states (DOS) con- 
tributions at the initial and final states momenta of the 
Brillouin zone which lie on a surface of constant energjPl 



, I.e. 



6N^{lo) (X / iV°(c.)iVO+q(c.)dk, 



(12) 



where N^{uj) is the spectral density function (DOS con- 
tribution) originating at momentum k at given bias en- 
ergy w. 

III. NUMERICAL RESULTS 

Now we start to examine the effect of the single im- 
purity scattering center on the quasiparticle interference 
pattern of CeCoIns. Since the corrugated Fermi surface 
of CeCoIns has considerable three dimensional character, 
we are not able to use the T-matrix formalism directly 
to compare with the two dimensional (surface) FT-STM 
results. In this respect we restrict ourself to the ab plane 
by averaging over th e momenta in the direction. 
Using the Eqs. ( [sfiol one can calculate the T-matrix 
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FIG. 3. (Color online) Total quasiparticle interference (quasi- 
particle interference in Fig. [2] averaged over fcz), for normal 
state (first panel) , superconducting state with d^^y gap sym- 
metry (second panel) and superconducting state with d^i^yi 
gap symmetry (third panel). (For bias voltage cj = O.OlAo 
and Ao/PF = 0.086). First row corresponds to the non- 
magnetic impurity, second row refers to magnetic impurity 
and last one shows the difference 5QPI between quasiparticle 
interference for magnetic and non-magnetic impurity scatter- 
ing. 
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FIG. 4. (Color online) Same as Fig. [s] but for bias voltage 
u = O.lAo. 
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and corresponding Green's function for each slice of the 
Fermi surface at a particular fc^, then by averaging over 
the individual QPI of each slice it is easy to find 
the total QPI for tunneling current along (OOl)-direction. 

The result of the QPI for the non-magnetic impu- 
rity is presented in Fig. [2j for individual layers at 
k^c = 0,0.27r,0.47r,0.67r,0.87r,7r. The first panel of this 
plot shows the corresponding quasiparticle interference 
strength in the normal state, and to compare the different 
superconducting gap symmetries in this compound, we 
look at the two main candidates discussed above, namely 
dxy gap symmetry, 

Ak = Aq sin k^a sin kyU, 

and dx-2_y-2 gap symmetry, 

Ak — Ao(cos kxO, — cos fc^a), 

in the second and third panels respectively. In Fig. [l|b) 
the constant energy intensity maps of the spectral density 
for normal state (first panel) superconducting state with 
dxy (second panel) and d3.2_j,2 (third panel) gap symme- 
try are plotted for the different cuts of the Fermi surface 
at kzC = 0, O.Stt, tt. Their node structure with respect to 
the Fermi surface is indicated by dashed lines. 

First we analyze the results for individual k^ cuts for 
nonmagnetic impurities. One can see that for the cases 
with kzC = 0, 0.27r, 0.47r, the gap with dxy symmetry does 
not have a node on the Fermi surface. Since the bias volt- 
age is smaller than the gap value, DOS is reduced and as 
a result the QPI becomes negligible. For larger k^ (be- 
cause of the inner structure of the spectral density func- 
tion) when the dxy gap becomes nodal or the node-lines 
approach the Fermi surface a different behavior of QPI 
appears. An enhancement of QPI for the distinguished 
q-vector which joins the new node points is observed for 
kzC = 0.67r,0.87r (see qi and q2 in Fig. [IJb.V)), by ap- 
proaching to k^c — t: the size of the inner part is growing 
thus the effect of the node and QPI become smaller again. 

This scenario is completely different for the dx2-y2 case 
where the node lines cut the Fermi surface at any value 
of kz- Because in the regions close to the node line, the 
DOS contributions of the scattered electrons are rapidly 
growing, QPI will be enhanced for the connecting q^ - 
vectors, and a corresponding point like pattern at q^ can 
be found in all values of kz- Especially when the inner 
part of spectral density function appears the number of 
these points, as shown in last panel of the Fig[T}D, are 
increased. 

The observable tunneling current is an average over 
the Fermi surface slices at different kz- Therefore we 
have to consider the calculated QPI averaged in (001)- 
direction. For this reason we plot in Fig. [3] the averaged 
QPI pattern, for normal state in the first panel, and su- 
perconducting state with dxy and dx^-y^ gap symmetry 
in second panel and third panel, respectively. 

In this figure the first row corresponds to QPI from 
non-magnetic impurity scattering which is calculated by 



averaging over the individual kz- layers presented in 
Fig. ([2]). In the normal state a trace of the dominant q 
- vectors which join the prominent points of the spectral 
density functions (Fig. [l]b first panel) can still be fol- 
lowed. The QPI is most pronounced at small q - vectors 
and also around the zone boundary vectors (±7r,±7r). 

In the superconducting state with dxy gap symmetry 
due to the gap opening the strength of the QPI decreases, 
and if we neglect the inner structure of spectral den- 
sity the remaining part has very small effect in QPI. By 
considering the inner part of the spectral density func- 
tion which has a node, pronounced QPI pattern can be 
found at corners (±7r, ±7r) and also in (100)- and (010)- 
directions. For a dx2_y2 gap symmetry QPI pattern in 
the Brillouin zone is distinctly different from normal state 
as well as dxy state. In this regime the pronounced points 
QPI pattern can be found at small q vectors in (110)- , 
(OlO)-and (lOO)-directions. 

The second row of Fig. [3] refers to QPI strength due to 
magnetic impurity scattering. The general behavior for 
the magnetic impurity is the same as the non magnetic 
impurity. Consequently we do not see dramatic change of 
QPI pattern aside from an overall reduction of the QPI 
amplitude at small momenta. Instead we focus on the 
most interesting QPI quantity for a comparison with pos- 
sible experimental results which is the difference of QPI 
strength for the magnetic (m) and non-magnetic (n.m) 
impurity scattering defined by i5QPI = QPIm — QPIn.m- 
We plot this difference in the last row of the Fig. |3] 
Magnetic impurity scattering leads to enhanced QPI for 
the large q- vectors, and to reduced QPI for the small 
q- vectors, in the superconducting states. The QPI in 
the superconductor with dxy gap symmetry is along the 
(lOO)-and (OlO)-directions and there are some effects at 
(±7r,±7r), but for the dx^^y-i gap symmetry we observe 
small q vector structures along the diagonal and we do 
not see any effect at large vectors around (±7r,±7r). 

By increasing the bias voltage size the small pockets 
in spectral density function are growing which enhances 
these effects. In the Fig. [4] we present the QPI for the 
larger bias voltage lo = O.IAq. As we expect by increas- 
ing the size of the electronic pockets, the QPI pattern 
grows, and the intensity shifts to higher momenta. The 
QPI still has an axis-aligned cross structure for dxy and 
diagonal cross continuous pattern with sharp points along 
(100)- and (OlO)-directions for c?^2_j^2. By increasing the 
bias voltage the QPI patterns which are created by the 
inner part of spectral density function are growing for 
both symmetries. However the sharp point like picture 
along (100)- and (OlO)-directions in for the dx^-y-i and 
at (±7r, ±7r) for dxy are still distinguishable. 

Finally we remark on the effect of a magnetic field act- 
ing only on the spin degrees of freedom. Our calculations 
show that adding the Zeeman term 7^5= u^bB to the 
quasiparticle dispersion does not effect the previous re- 
sult dramatically for accessible field strength. Another 
more promising direction would be to study the influ- 
ence of a FFLO type spatially inhomogeneous supercon- 
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ducting order parameter on QPI pattern in the high field 
regime where the FFLO state has been identified e.g. in 
recent NMR experiments^. 

IV. SUMMARY 

In this paper we presented the theory for the SI-STM 
in the heavy Fermion superconductor. We study the 
quasiparticle interference in heavy Fermion superconduc- 
tor CeCoIns systematically. We compare the effect of the 
different singlet superconducting order parameters that 
have been proposed, namely dxy and d.j;2_y2 , on the QPI 
in this compound. By comparing the QPI pattern of sin- 
gle non-magnetic and magnetic impurities we have shown 
that QPI has an axis aligned cross structure with some 
features at (±7r,±7r) for dxy gap function, and a diago- 



nal cross continuos pattern with sharp points along (100)- 
and (OlO)-directions for dx2-y2 gap symmetry. Further- 
more we have shown the difference of QPI in present 
of magnetic and non-magnetic impurity, (5QPI, becomes 
positive for large and negative for small wave vectors. 
We have shown that the axis aligned and diagonal cross 
structures are clearly seen in the (5QPI for dx2_y-2 and 
dxy gap symmetries respectively. 

We conclude that QPI pattern are sensitive to the or- 
der parameter symmetry in the 115 compounds and may 
be observed in future STM measurements. This obser- 
vation would provide a direct fingerprint to distinguish 
between different superconducting order parameter mod- 
els. It may be particularly relevant for the case of Celrlns 
where the symmetry of the supe rconducting gap function 
is discussed controversially^SEZl, 
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